Are Directed Waves Multifractal? 
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Wave propagation is studied in a sufficiently anisotropic random medium that backscattering 
along one direction can be neglected. A Fokker-Planck equation is derived the solution to which 
would provide a complete statistical description of such directed waves. The Fokker-Planck equation 
is mapped onto an su(l,l) ferromagnet and its symmetries are identified. Using the symmetries 
asymptotic wave function distributions are computed and used to show that directed wave functions 
fill space uniformly and do not have multifractal character. 
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Much progress has been achieved in the study of wave 
propagation in random media by the infusion of scal- 
ing ideas from statistical mechanics and critical phenom- 
ena. Dirty electronic systems are the best studied ex- 
ample 0,^]. They are now known to exhibit a variety 
of phases and transitions, many of which remain poorly 
understood. Recently there has been interest in the ap- 
plication of multifractal analysis to these problems 
Fractals are complex geometric forms that exhibit self- 
similarity under magnification || . Multifractals are still 
more complex objects which need an infinite set of expo- 
nents to characterise their self-similar scaling rather than 
the single exponent that suffices for an ordinary fractal 
j|. These ideas will be expressed precisely below; no 
previous background in multifractal analysis is needed to 
read this paper. 

In statistical mechanics the typical configuration of a 
disordered spin-system at a critical point is complicated 
and for certain models can be shown to exhibit multi- 
fractal scaling The purpose of such an analysis is to 
provide a characterisation of the critical point in terms of 
multifractal exponents that supplement the conventional 
universal exponents of critical phenomena Q . Dirty elec- 
tronic wave functions have also been studied from this 
perspective. For example, extensive numerical simula- 
tions indicate that the wave functions of a quantum Hall 
system are multifractal at the critical point that sep- 
arates two Hall resistance plateaus Calculation of 
multifractal exponents requires detailed statistical infor- 
mation about the wave functions and analytic results are 
therefore difficult to obtain. With few exceptions 0, 
multifractal analysis for dirty electronic systems is there- 
fore carried out via numerical simulation. In this Letter 
we report multifractal analysis of directed wave propaga- 
tion through a random medium. This problem has been 
investigated by the authors of ref |8) who claimed numer- 
ical evidence of multifractal scaling. 

The directed wave model is very rich in physical ap- 
plications. It applies whenever waves propagate in a suf- 
ficiently anisotropic medium that backscattering in one 



direction can be neglected PjuJ. An experimental real- 
ization |Tl|-[l^| that has attracted much attention recently 
Jld] po[ (and refs therein) is the surface of a quantum Hall 
multilayer. Studies of directed waves have a long history 
due to the equivalence of this problem to time-dependent 
wave propagation in a noisy environment pl|j8|. More- 
over a recent model for stress propagation in granular 
materials p2] is identical to lattice models of directed 
waves p3| . Still another motivation for studying the di- 
rected wave model has been its formal similarity to the 
directed polymer problem pij (for an illuminating dis- 
cussion of this viewpoint see, for example, ref |2j|). 

Previous analytic work on directed waves has focussed 
on calculation of low-order correlations of the wave func- 
tions (needed, for example, to calculate observables such 
as the disorder averaged conductance of a quantum Hall 
multilayer) . Here it is neccessary to analyse the full prob- 
ability distribution of the wavefunctions. Our method is 
to derive a Fokker-Planck equation that governs the evo- 
lution of the probability distribution. This equation has 
an su(l,l) symmetry that can be exploited to obtain the 
asymptotic probability distribution which reveals that di- 
rected waves fill space quite uniformly and therefore do 
not change appearance under magnification. Thus, con- 
trary to previous numerical work, we are able to show 
exactly that directed wavefunctions do not have a mul- 
tifractal character. The asymptotic wave function distri- 
bution (eq 3 below) is the directed wave analogue of a 
result derived several years ago for directed polymers in 
one-dimension by Huse et al. |27| . 

We shall take directed waves at a fixed frequency to 
evolve according to (see for example, [ p7| , ^8| ) 

d 

- l ~Q^^ n ( X ) +tn(x)4>n+ l(x) +t n _ 1 (x)'il) n -l(x) = 0. 

(1) 

Backscattering is neglected along the x— direction. The 
transverse direction is taken to be discrete and the system 
is assumed to be of finite size, N, in this direction (with 
periodic boundary conditions). The anisotropic and di- 
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rectional nature of eq (1) is reflected in the fact that it is 
of first order in the x— direction whereas it is of second 
order in the transverse direction. Disorder is incorpo- 
rated by taking the hopping terms t n (x) to be random 
with statistics given below. Other models have been con- 
sidered in the literature. They are either equivalent to 
eq (1) or are believed to have the same qualitative be- 
haviour [^6|. The problem posed is the following: The 
wavefunction is specified for a fixed value of x; it is then 
evolved in the x— direction according to eq (1). For def- 
initeness, it may be supposed that the wavefunction is 
localized at n — at x = 0. Thus ip n (x = 0) = <5„.o- 
Each disorder realization will produce a different wave 
function at a larger value of x > 0. We are interested in 
P{a n ,x), the probability density that at a fixed x > 
the wave function is ipn(x) — d n . 

Our problem resembles the problem of Brownian mo- 
tion, in which a heavily damped particle is subject to a 
noisy environment. Standard methods exist in the theory 
of Brownian motion to convert the stochastic equation 
of motion (called the Langevin equation) into the corre- 
sponding Fokker-Planck equation which governs the time 
evolution of the probability density of the particle posi- 
tion |p8| . These methods may be applied to obtain the 
Fokker-Planck equation corresponding to eq (1) 

-dP{a n ,x)/dx = HP; 
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where the Greek letters represent either 1 or 2. Here 
we have refined our notation to write a 1 = real part 
of the wave function and a 2 — imaginary part. Like- 
wise t 1 and t 2 denote the real and imaginary parts of 
the hopping term. In deriving eq (2) it was assumed 
that the hopping is a Gaussian white noise process with 
zero mean and variance given by [i"(x)i^(a;')] ons = 
{D /2)5 a ,f)5n,m5{x — x'). [. . .] ons denotes an average over 
the ensemble of disorder realizations. Given an initial 
wavefunction, ip n (x = 0), the Fokker Planck eq (2) in 
principle allows calculation of the complete probability 
distribution of the wave function, P(a,x) for larger val- 
ues of x. In practice, this may appear difficult because 
eq (2) is a partial differential equation in a large number 
(2N+1) of variables. However it will be seen that eq (2) 
has a number of helpful symmetries. 

First note that stationary solutions to the Fokker 
Planck equation must be radial. It is easy to verify 



by substitution that if P is radial (a function only of 
J2 n a ( a n) 2 an d possibly x) then dP/dx = showing that 
radial solutions are stationary. The converse, that sta- 
tionary solutions must be radial, is also true and will be 
shown below. It is a consequence of an su(l,l) symmetry 
of the Fokker-Planck equation to which we shall return. 

A second useful property of eq (2) follows from prob- 
ability conservation. By direct substitution into the 
wave eq (1), it can be shown that the total probabil- 
ity J2 n a(V'ra) 2 does not change with x for any disorder 
realization. Consequently if we solve the Fokker-Planck 
equation subject to the initial condition that the system 
begins with some definite normalized wave function, P 
must live on the unit sphere in a— space. In other words, 
P vanishes unless J2 n a {a") 2 = 1. 

The only stationary distribution consistent with prob- 
ability conservation is a uniform distribution on the unit 
sphere in a— space. Thus the large x asymptotic distri- 
bution is uniquely determined to be 
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Here A is a constant fixed by the normalization condition 
I Yin a da^P(a, oo) = 1. The physical content of eq (3) 
is that directed wave functions are completely random- 
ized after propagating a large distance; however they re- 
main normalized. Precisely the same distribution arises 
in connection with the distribution of matrix elements 
of correlated random matrices |p9f ; thus techniques for 
averaging over such distributions are well developed. 

Multifractal characterization of the directed wave func- 
tions begins with calculation of the exponents C(<?) which 
are defined via B 
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The generalized fractal dimensions D(q) are related to 
these exponents by ((q) = (1 — q)D(q). For q — 2, the 
left hand side of eq (4) is the inverse of the participation 
ratio which has long been studied in connection with dis- 
ordered electronic systems [M . Roughly it measures the 
number of sites on which the wave function has a sub- 
stantial weight. Eq (4) shows that £(#) characterizes the 
growth of generalized inverse participation ratios with 
system size. Evidently £(0) = 1 and £(1) = (by nor- 
malization). To calculate the multifractal dimensions for 
other values of q, it is neccessary to average the inverse 
participation ratios over the wave function distribution 
in eq (3). 

Evaluation of the necessary integrals is facilitated by a 
trick given in ref p9| . For illustration consider the nor- 
malization integral A — J daS(a 2 — 1). Here for brevity 

Iln a d a n ^ s written as da; J2 n a ( a n) 2 as fl2 ; an d in- 
dices are suppressed everywhere. Rescale by introducing 
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b = ay/r and find r N ~ l A = J dM(6 2 - r). If both sides 
are multiplied by e~ r and integrated over the rescaling 
factor r from zero to infinity, the remaining b integrals 
become straightforward gaussians and the final result is 
A = tt n /T{N). The other integrals can be evaluated in 
a similar fashion to yield 
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(5) 



The last approximate equality holds for large N. Com- 
paring eq (4) and (5) we see that = 1 — q (which 
implies D(q) = 1). This is precisely the scaling expected 
of an object that is not multifractal. 

To make contact with previous numerical work [pi it 
is neccessary to calculate the g((3) spectrum which is 
defined as follows Let w be the weight (modu- 

lus square of the wave function) on a particular site. 
f3 is a logarithmic measure of the weight defined as 
(3 = In w/ In TV. Assemble a histogram of (3 values by 
drawing from the ensemble of disorder realisations and 
let HN(f3')d[3 represent the frequency with which (3 lies 
between f3' and j3'+df3 in a system with N sites. For large 
N it is expected that H N ((3) ~ N 9 ^ which defines g{/3). 
To calculate g(/3) it is helpful to first evaluate P(w)dw, 
the probability that the weight on a particular site lies 
between w and w + dw, obtained by integrating the dis- 
tribution in eq (3) over all sites except one. For large 
N the result is P(w) = N exp(-Nw). Straightforward 
substitution into the definition above then yields 



?(/?) = 2 + /3 



(6) 



for (3 < -1. Note g(0) > over the range -2 < (3 < -1. 

In ref || the g((3) spectrum is calculated numerically 
and found to have support on the interval —1.97 < (3 < 
—0.88. The authors of ref || assert (incorrectly) that if 
g(/3) has support over a range away from (3 — —1, it nec- 
cessarily implies multifractal scaling. This is the primary 
basis for their claim of multifractal scaling as the actual 
deviation from ordinary scaling they infer is very small: 
((q) (numerical) = 1 — q + 3.6 x 10 — 3 g 2 . The asymptotic 
distribution derived here, eq (3), explicitly shows both 
ordinary scaling and a spectrum, g((3), with an extended 
support away from [3 = — 1. Moreover the form of g{(3) 
we obtain, eq (6), provides a good fit to the numerical 
data (fig 3 of ref ||). Hence we believe that the numerical 
data of ref j^] are in fact consistent with our conclusion 
that directed waves are not multifractal. 

To complete the analysis we must now return to the 
Fokker-Planck equation and show that all stationary so- 
lutions are radial. Note that eq (2) has the appearance 
of an imaginary time Schrodinger equation. From this 
point of view P(a, x) is the wave function and we are in- 
terested in the multiplet of eigenfunctions of the "Hamil- 
tonian", Ti., with eigenvalue zero. It is useful to express 



7i in second quantized language by introducing harmonic 
oscillator ladder operators 

b? = 4= f < + IT-) ; (KV = "4= (< - 77~ — J ; (7) 



which obey the bosonic commutation relations [6", i 
SafiSnm and [6",6^J = 0. In this language the Hamilto- 
nian is given by 
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Thus the Fokker-Planck equation is seen to be equivalent 
to an interacting boson problem. 

The symmetry of the Hamiltonian is revealed by con- 
sidering the algebra of the bosonic bilinear terms out 
of which it is built. Define K+ = i^i^j', K~ = 

\ E a KK, K* n = \{b^bl + bl&J) and M n = i^Jbl - 
b 2 n b\). These operators obey 



[K+,K-] = -2K z ; \K\K ± \=±K i 



(9) 



and M commutes with all the ICs. In writing eq (9), 
the site indices have been suppressed for clarity. These 
commutation relations apply only if the indices coincide. 
Operators with distinct indices all commute. This alge- 
bra is reminiscent of the angular momentum algebra— the 
notation was chosen to highlight the similarity. However 
there is a crucial sign difference in the [K + , K~\ commu- 
tator. Thus the K operators actually obey the su(l,l) 
algebra rather than su(2). su(l,l) is a well studied clas- 
sical Lie algebra sometimes called the hyperbolic angular 
momentum algebra in the physics literature when it is 
discussed in connection with Schwinger's coupled-boson 
description of angular momentum pl| ]. 

The Hamiltonian may be rewritten in terms of these 
operators as an su(l,l) ferromagnet (a generalization of 
the ordinary Heisenberg ferromagnet in which su(l,l) op- 
erators replace their spin counterparts) [p2| : 



H = DJ2 [2KK+1 - (K+K-+1 + K-K+ +l ) 

n 

-l(M n M n+1 + l)] 



(10) 



It is now easy to verify that the Hamiltonian is invariant 
under su(l,l) rotations: 



[H, Kt, 







(11) 



3 



where K tot = J2n K n 

The complete set of eigenstates with eigenvalue zero 
may now be found. On account of the symmetry ex- 
pressed in eq (11), this multiplet should form an irre- 
ducible representation of the su(l,l) algebra. Note that 
the boson vacuum is a zero energy eigenstate of the 
Hamiltonian. This is evident from inspection of eq (8). 
In addition, the vacuum is annihilated by K^ ot . Thus 
the vacuum is the lowest weight state in the multiplet; 
the complete infinite set may be obtained by repeatedly 
applying K* ot to the vacuum pl[ . 

To see that each member of the multiplet (K^ ot ) l \0 > 
with I — 0,1,2,... corresponds to a radial function it 
is neccessary to translate back into first quantized lan- 
guage. The ground state of a harmonic oscillator has a 
gaussian wave function; hence the boson vacuum is given 
by 



|0>-ex P -i]T«) 2 



(12) 



-a radial function. K^ ot is given by 
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It is easy to verify that the action of K^ ot on any radial 
function will yield another radial function. Thus we have 
shown that the vacuum is radial as is any state built out 
of it by repeated application of the raising operator K+ ot \ 
in other words, a stationary solution must be radial. 

Finally, note that although in this paper we have fo- 
cussed on stationary solutions, it may be possible to give 
a rather complete analysis of the directed wave problem 
due to the equivalence demonstrated here to a relatively 
simple model— an su(l,l) ferromagnet. 

In summary the main results of this paper are: deriva- 
tion of a Fokker-Planck equation that governs the wave 
function distribution for directed waves; mapping of this 
equation onto an su(l,l) spin-chain which reveals its sym- 
metries; and calculation of the asymptotic wave function 
distribution which shows that directed waves fill space 
quite uniformly and are not multifractal. 
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